Memory effects in non-interacting mesoscopic transport 
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Abstract 

Consider a quantum dot coupled to two semi-infinite one-dimensional leads at thermal equilibrium. 
We turn on adiabatically a bias between the leads such that there exists exactly one discrete eigenvalue 
both at the beginning and at the end of the switching procedure. It is shown that the expectation on 
the final bound state strongly depends on the history of the switching procedure. On the contrary, 
the contribution to the final steady-state corresponding to the continuous spectrum has no memory, 
and only depends on the initial and final values of the bias. 

1 Introduction and the main results 

Memory effects in quantum transport are rather common in systems where the carriers have self- 
interactions, and they can be of several types. Both in the Master Equation approach (T5l [16] and 
the TDFT approach [22] one observes a dependence of the steady state current both on the initial state 
and the switching procedure. Very recently [5] such a dependence on the initial state of the sample was 
proved in the cotunneling regime. 

In contrast, in the non- interacting mesoscopic quantum transport one can show that the charge/energy 
current observables depend neither on the initial state of the sample, nor on the switching procedure. 

The sudden switch of a coupling at t = (in an initially partitioned system) has been thoroughly 
investigated in a number of previous works, see for example [U HU Q21 [l3j [19]; for a more physical 
approach see [3J. The steady state currents (computed as Cesaro means) are very robust, irrespective of 
which method one uses in order to induce a non-equilibrium state in the system (see |10j): their values 
only depend on the initial equilibrium state and on the final expression of the Hamiltonian which governs 
the evolution after the switching procedure is over. 

In the partition- free approach introduced by Cini [4], the situation is rather similar. In [5] it is 
treated a situation in which a thermal equilibrium state is perturbed by turning on adiabatically a bias 
between the leads. It is shown that if the instantaneous discrete spectrum of the one-body Hamiltonian 
is always well separated from the instantaneous continuous spectrum, then the adiabatic limit of the 
current coincides with the steady state current value of the sudden switch f [2"6l |2"?1 lo] ) . 

In the current paper we investigate the situation in which the instantaneous discrete spectrum can 
enter the continuous one. It turns out that at the adiabatic limit, the steady state value of the charge 
current still has no memory of the switching procedure. But the situation is totally different for the 
expectation on the final bound states. The adiabatic limit of this expectation is highly dependent on 
whether the instantaneous discrete spectrum enters the continuous spectrum or not. From a mathematical 
point of view, this phenomenon is related with the so-called 'adiabatic pair creation', see [20 ] |2T ] 124]. 

Our model is of Wigner-Weisskopf type, see [H] for a rather complete spectral analysis. We choose 
this model because we want to maximize clarity and minimize the technicality of the proofs, but most of 
the results below can be generalized to samples containing more sites, or even to a continuous setting. 

1.1 The setting and notation 

We consider two semi-infinite discrete leads coupled to a small system consisting of just one site. 

The single-particle Hilbert space is H — C © {Z 2 (N_) © Z 2 (N+)} =: Hs © Hl- The canonical basis 
in Hl is denoted by {|i 7 ) : 7 = ±, i > 0} where i 1 is the i-th site of the lead 7. Similarly, we denote 
by {\S}} the basis element of C. With these notations we introduce the single-particle Hamiltonians 
h± which describe an electron on the leads to be just two copies of the usual one-dimensional discrete 
Laplacean L initially defined on l 2 (Z) and then restricted to Z 2 (N) with Dirichlet boundary condition at 
— 1. With a physicist's notation we have 

L = + l)<j| + l)<j|}. 
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If Il-t = X^>o b±)0±l are the projections on the left/right leads, then by definition h± := II-1-LII-1-. We 
introduce the operators: 

h L :=Y,K> h s :=E \S)(S\, h T = r £{|0 7 )(S| + |S)(0 7 |}, (1.1) 

7— ± 7— ± 

where Eq > and < |t| < 1 are real parameters to be chosen later. If v > is another real parameter 
which models the potential bias between the two leads, then the total Hamiltonian reads as: 

h(v) :=h s + h L +vn_ +h T . (1.2) 

We write 

ho(v) = h s + h L +vIl-, (1.3) 

such that h(v) = ho(v) + Iit- 

The spectrum of h± is absolutely continuous and equals [—2,2]. The continuous spectrum of ho(v) 
is <y ac (h (v)) = [—2,2] U [—2 + v,2 + v] while the pure point part is independent of v and given by 

fppCioO)) = {^o}- 

In Section[2]we will treat in great detail the spectral properties of h{v) as a function of v. In particular, 
we will show that if Eg is large enough and r small enough, then h(v) has a unique discrete eigenvalue 
X(v) in the interval (v + 2, oo) as long as v is strictly smaller than a critical value w Cj i. When v — v Cj \, 



there is no point spectrum at the threshold v c ,i + 2 (see Proposition 2.1 ). Moreover, there exists a second 
critical value v c 2 > w c ,i such that if v G [v c ,i, v c ,2] the spectrum of h(v) is purely absolutely continuous, 
but if v > v Ct 2, an eigenvalue X(v) appears in the interval (2, v — 2) and stays there. 

In order to model the switching procedure, we need to make v time dependent. We will only consider 
switching procedures for which v(s) — if s < — 1. Moreover, we make the following assumption: 

Assumption, v is twice piecewise differentiable on (—1,0) with uniformly bounded second derivative, 
and v is continuous at —1 and 0. 

We are interested in situations where the discrete spectra of h(v(—l)) and h(v(0)) consist of precisely 
one eigenvalue. We will only consider the following three generic situations. 

1. The first case is when the potential bias v is C 2 on [—1,0] and does not cross the critical values, 
such that the discrete instantaneous eigenvalue is always present. See Figure [l] 

2. The second situation is when the bias potential crosses twice the critical value v c i, causing the 
instantaneous eigenvalue X(v(s)) to disappear at some point in time and then to reappear at a later 
moment. We model this potential bias by a function v : [—1,0] — > [0,u c ,2)) continuous at —1 and 0, 
with u(— 1) = and v(Q) = w C) i — 1. We assume that there exist — 1 < s c < s' c < and < 6 << 1 
such that: 

v(s c - 0) = v(s' c + 0) = v Ct i - S, 

v(s c ) = v(s c + 0) = v c>1 +6 = v{s' c ) = v(s' c - 0), 
v c ,i + S < v(s) < v Ci2 , s G [s c , s'J; 

v(s) < v c<1 -5, s<£ [s c ,s' c ]. (1.4) 

See Figure [2] 

Keeping in mind that the coupling constant r in Kt must be small, this potential bias insures that the 
instantaneous Hamiltonian h(v(s)) will have exactly one discrete eigenvalue if s G [—1, s c ) U (s' c , 0], 
and purely absolutely continuous spectrum for s G [s c , s' c ]. At the end we will let 5 go to zero, but 
until then the potential bias has a small discontinuous jump at s c and s' c . 

3. The third physically interesting switching procedure (from the mathematical point of view being 
though closely related to the second situation) is the one in which the bias causes the instantaneous 
eigenvalue X(v) to disappear into the continuous band [—2 + v, 2 + v] for u c l < v < w c 2 and to 
reappear and stay in the interval (2,v — 2) for v > v Ct 2- More precisely, we consider an increasing 
function v : [—1, 0] — > [0, v C: 2 + 1], continuous at —1 and with v(— 1) = and v(0) = w c ,2 + 1< We 
assume that there exist — 1 < s c < s' c < and < 6 << 1 such that: 

v(s c - 0) = v c> \ - <5, v(s' c + 0) = v c . 2 + 5, 

v(s c ) = v(s c + 0) = v c ,i + S, v(s' c ) = v(s' c - 0) = V c ,2 - S. (1.5) 

See Figure [3] 
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Figure 1: The first situation 
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Figure 2: The second situation 




Figure 3: The third situation 
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In what follows, we will replace the notation h(v(s)) by h(s) at all points where v is continuous, and 
also write h(s ± 0) := h(v(s ± 0)) for the (finite) set of points s at which v has discontinuities. Moreover, 
we will adopt the same convention for various functions of h(s) e.g.: 

Pd(s) = P d (h(s)) = \i>{a))(il>(s)\, P ac (s) = P ac (h(s)), (1.6) 

where ip(s) is an eigenstate corresponding to the (only one) eigenvalues of h(s). We recall that P c(^l) = 
II_ + n + . Let us point out that the wave operators 

n(s) := s-lim e ith ^e- ith ^P ac (h L ) (1.7) 

t— > — oo 

exist and are asymptotically complete [25] , 
1.2 The results 

If rj > is a small adiabatic parameter, we consider the time dependent Hamiltonian 

h{v[rjt)) =h s + h L + v(r)t)H- + h T = h (v(r]t)) + h T . 

We denote by U(t,to) the unitary solution of the time dependent Schrodinger equation 

iU'(t,t ) = h(r)t)U(t,t ), U{t Q ,t Q ) = l, -rf x < t, t < 0. 

The initial equilibrium state of the system is characterized by a density matrix operator p eq which is 
assumed to be a function of the operator h(— 1) = ho + hr, i.e. p cq = f eq (h(— 1)). An example could be 
/eq(x) = l/{e^ x -^ + 1) with /3 > and /x € R. For physical reasons we demand that < ||/ cq ||oo < 1< 
Note that the one-particle density matrix is not a trace class operator. Given a trace class observable A, 
its equilibrium expectation is given by Tr{p cq A}. 

In the remote past t < —r/^ 1 , the density matrix operator equals p eq and it is time independent. At 
time — we start out the potential bias v and we let it evolve. 

The density matrix operator solves the Liouville equation ipL(t) = [h(r]t), p v (t)] on the interval —rj" 1 < 
t < 0, and it is given by the formula 

p n (t) ^U^-rj-^p^U^-r,- 1 )*. 

At t = 0, our state will be described by: 

/>„ := Pv (0) = C/(0, -^WE^O, -T 1 )* = U(- V -\Q)* p cq U(- V -\Q), 

where in the last equality we used the fact that U(t, t')* — U(t',i). 

If A 6 P>i(H) is a given self-adjoint trace class observable, the question we want to answer in all three 
cases is the existence of the following adiabatic limit: 

(A) := lim Tr{p n A} = lim Tr{U{-r]-\ 0)* f^hpi-rj- 1 , 0)A}. (1.8) 

If A = P<j(0) = |^(0))(V'(0)|, then (^4) represents the probability of arriving at the final bound state at 
the end of the switching procedure. Now here is our main result: 

Theorem 1.1. We have the following situations: 

(i) If the potential bias v is chosen in such a way that the instantaneous eigenvalue X(s) of h(s) is always 
bounded away from the instantaneous continuous spectrum and stays in the interval (v(s) + 2, oo), then: 

(A) = Ti{n(o)f cq (h L )p ac (h L )n(oyA} + (?/>(o) | a^(o)) / eq (A(-i)). (1.9) 

If A = P d (0) = |V(0))<V(0)|, then (P d (0)) = / oq (A(-l)). 

(ii) If the potential bias v satisfies either (1.4) or (1.5) then 

(A) = Tr{fi(0)/ oq (/i L )P ac (/i L )r!(0)M} + (V>(0) | At/>(0)) 

■ Ms' c + o) | n(s' c - o)f cq (h L )p ac (h L )n( s ' c - 0)^(4 + o)). (1.10) 

Moreover, 

(p d (o)> - + o) | n(s' c - o)f cq (h L )p, c (h L )n( s ' c - o)>(4 + o)). (1.11) 
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Remark 1.2. The first terms in both ( 1.9 ) and ( 1.10 ) are identical and only depend upon the initial and 
final values of the bias potential. Moreover, they are the only ones contributing to the adiabatic charge 
current since in that case A = i[h(v(0)),H±] which implies (V'(O) | ^4-^(0)) = 0, thus the second terms 
disappear. This is consistent with previous results [10] showing that the adiabatic limit of the current 
does not depend upon the specific form of the switching. On the contrary, while in (i) the second term 
in the r.h.s. of ( 1.9 )) again only depends upon the initial and final values of the bias potential, in (ii) the 
second term in the r.h.s. of ( 1.10 ) depends upon the switching procedure but only via the behavior of the 
bias potential in the neighborhood of s' c , i.e. at the last passage through a critical value. The fomulas 



(1.9 1 and (1.9) are explicit, but the price is that we assumed a jump of size 26 for the bias potential at 



the critical values. 

Our second main result gives (Pd(0)) for all three situations in the limit 5 \ 0, i.e. 
discontinuity jump in the bias potential shrinks to zero. 



when the 



Proposition 1.3. (i) If the bias v varies in such a way that the discrete eigenvalue X(s) of h(s) always 
remains separated from the continuous spectrum (see Figure^, then: 



<Pd(0)) = / eq (A(-l)). 



(1.12) 



(ii) If X(s) enters the continuous spectrum but reappears and stays in the interval (v(s) + 2,oo) (see Figure 
[1]), we have: 

lim(P d (0)) = / cq (2). (1.13) 

(ii) If X(s) enters the continuous spectrum but reappears and stays in the interval (2,v(s) — 2) (see Figure 
[M, we have: 



lim<P d (0)) = / cq (-2). 



(1.14) 



Remark 1.4. It might happen that / oq (A(— 1)) = 0, i.e. the initial state has no discrete component. If 
the continuous spectrum is not crossed, the probability of finding the system in the final discrete state 
is still zero (see (1.12)). But if either / eq (2) = 1 or / cq (— 2) = 1, by entering the continuous band this 



probability can be made equal to 1. 

One might think that Proposition 1.3 (ii) and (iii) also cover the generic (i.e. without jumps) bias 
potentials. Unfortunately this is not true since we first took the adiabatic limit rj \, 0, and after that 
S \ 0. In order to cover the general case one has to perform the limits in the reversed order. As it is 
known from the spontaneous pair creation case PHI 1211 121] this is a hard technical problem (even more 
demanding for the problem at hand) since one has to control the evolution in the adiabatic limit near the 
critical times and this relies on detailed (model dependent) spectral and scattering study near criticality. 
Notice that contrary to the spontaneous pair creation problem, in our case one has to deal with a critical 
Hamiltonian having a resonance at the threshold, and not a bound state. We believe that a careful 
study of spectral and scattering theory near energetic thresholds enlarging and streamlining the results 
in [231 [Ml US US] is needed in order to solve this very interesting open problem. 

We can allow multiple crossings of both critical values, and the results are rather similar. The 
probability of finding the system on the final discrete state after taking the limit S \ will be given by 
/ eq (±2), depending on which critical value was the last one crossed. 



1.3 The contents of the paper 

After this introductory section which also included the main results of our paper, we continue in Section 
[2] with a detailed spectral analysis of the operator h(v). The results are rather straigthforward but are 
needed in order to study the behavior of the discrete eigenvalue around the critical values of the bias. In 



Propositions 2.1 and 2.2 we show that the eigenvector corresponding to this discrete eigenvalue becomes 
more and more delocalized in the leads as we get closer and closer to the critical values. 

In Section [3] we prove Theorem [TTTj The main technical tool can be found in Proposition [3TT] which in 
some sense can be considered to be an adiabatic theorem for the continuous spectrum (see [21 HH1 [IZl [7] 
for results about adiabatic limit in scattering theory). 

In Section[4]we prove Proposition 1.3 More precisely, we investigate the limit S \ of the probability 
of finding the system in the final bound state. The fact that the discrete eigenvector becomes very 
delocalized when we approach the critical values of the bias plays a crucial role in the proof. 
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2 Spectral analysis of h(v) 



Introduce the notation R(z) for the inverse of h(v) — z in H, and Rl{z) for the inverse of + — z 
in extended by on Hs, i-e. 

Rl{z) = 0®(h_+v- z)' 1 © (h+ - z)- 1 . 

Note that the v dependence is not made explicit in this notation. 
Let us split the tunneling Hamiltonian hx as follows: 



hr = h 



LS T nSLi 



h LS = T\Q_)(S\+T\0 + )(S\=h 



SL- 



(2.1) 



Define the effective Hamiltonian as h e g(z) := (Eq — z)\S)(S\ — /isl-Rl(z)/ils; when restricted to T-Ls, the 
effective Hamiltonian is h e g(z) = G(z;v)\S)(S\, where 



G(z;v) :=E -z- r 2 (0_ \(h-+v- z) _1 0_) - t 2 (0+ | (h+ - z) _1 + }. 
Then the Feshbach formula for R(z) reads (see e.g [8]): 

1 



R(z)=R L (z) + 



-{1 - R L (z)h LS }\S)(S\{l - h SL R L (z)}. 



(2.2) 



(2.3) 



G(z;v) 

Let us denote by r±(z) = (h± — z)^ 1 the resolvents of the leads extended by outside ^ 2 (N±) in the 
obvious way. With this notation (2.3| writes: 

1 



R{z)=R L (z) + 



\S)(S\ 



G{z;v)> /x 1 G{z;v) 
: \S){0+\r+(z)- 



\S){0.\r.(z-v) 



-r-(z-v)\0-)(S\ 

T 



G(z;v) G(z;v 
+ r_(z- v) |0_) (0_ |r_ (z - v) + r_ (z - v) |0_) (0 + |r+ (z) 

— ^- [- -r+ (z) |0+) (5| + r+ (z) |0+) (0_ |r_ (z - v) 
+ r+(z)|0 + )(0+|r + (z) 



What is needed in (2.2 1 and ( |2.4[ ) are the formulas for r±(z). We recall them from [5]. 
Let 

Ci(z) = |(l-0^47^). 



(2.4) 



(2.5) 



where \Jx is chosen to have the cut (— oo, 0) i.e. \/x := y'jxje^/ 2 ) ar s( x ) j arg(x) e (— 7r,7r). Notice that 
Ci(z) is real analytic in C\ [—2, 2] and vanishes at infinity like 1/z. Also, £i is decreasing on the separate 
intervals (— oo, — 2) and (2, oo), while Ci(±2) = ±1. 
As it was shown in [5] we have: 



(0±,(h±-z)- 1 ± ) =-&(«), 
Ci(*) 



(m±,(/i±-z) n±) = 



Ci 2 (^) 



I (CiW |m -" l -Ci(2) m+n+2 ), n,m>0. 



In particular (2.2) becomes: 



G(«;v) =-B -^ + r 2 {Ci(z-w)+CiW}- 



(2.6) 



(2.7) 



2.1 The point spectrum of h(v) for v > 

As it is well known from the Feshbach lemma, the discrete spectrum of the coupled operator h(v) is given 
by the real solutions of the equation G(x; v) — 0, where x [—2, 2] U [—2 + v, 2 + v}. In the introduction 
we announced that we are interested in the case when Eq is large and far from the unbiased continuous 
spectrum, while the coupling t is weak. Let us assume that Eq > 10 and < |r| << 1. 
For any fixed v > 0, by differentiating in (2.2) we obtain: 



d x G(x;v) < -1, x £ [-2,2] U [-2 + v,2 + v]. 



(2.8) 



G 



which means that the map G(-;v) is strictly decreasing on the intervals (—00, —2) and (2 + v,oo). If 
v > 4 then G(-; v) is also decreasing on (2, —2 + v). In addition, 



lim G(x;v) =T<x>. (2.9) 

:— >-±oo 



We first show that there are no eigenvalues in the interval (—00, —2). Indeed, using (2.9 1 and the inequality 
G{— 2; v) > E Q + 2 -t 2 sup l/>0 |Ci(— 2 — v) +Ci( — 2)| > (remember that r is small enough and E > 10), 
we conclude that G does not change sign on (—00, —2), hence there are no discrete eigenvalues there. 
Now consider the interval (2 + v, 00). Because G(-, v) is strictly decreasing, we have that: 

-00 < G(x; v) < G(2 +v;v), 2 + v < x < 00. 



For all v > 0, formula (2.7) gives: 

G(2 + v;v) ^ E - 2 - v + T 2 {d(2) + Ci(v + 2)}. (2.10) 

There exists a (unique) eigenvalue if and only if G(2 + v; v) > 0, thus we need to investigate G(2 + v; v) 
as a function of v. 

Because £1 is decreasing on (2, 00) we have that G(2 + v; v) is strictly decreasing with v. If < v < 4 
then G(2 + v; v) > 4 - t 2 sup [/ >Q |Ci(2) + C,i(y + 2)| > for small enough t. It means that there exists 
exactly one non-degenerate eigenvalue X(v) £ (v + 2, 00) if < v < 4. If u increases even more, then 
G(2+i>; u) decreases approaching zero. There will be a critical value i> Cj i > 4 such that G(2+v Cy i; v C; i) = 
and G(2 + v; v) < if v > v Ct \. The value of v Ct i must be close to Eq — 2 if r is small. We conclude that 
if v > v C) i there is no discrete spectrum on (v + 2, 00). 

Third, we need to investigate what happens in the interval (2, —2 + v) when v > 4. We see that 
G(2; v) > for all v > 4, and we have the inequality: 

G(-2 + v; v) < G(x; v) < G(2; v), 2 < x < -2 + v. 

Thus we only need to investigate the sign of 

G(-2 + v; v) = E + 2 - v + t 2 {Ci(-2) + Ci(-2 + v)}, v > 4. 

Note that G(— 2+v; v) strictly decreases with v. If r is small enough, then G(—2+v; v) > G(2+v; v) > for 
all 4 < u < u Cj i, which means that there is no discrete spectrum in the interval (2, — 2 + v) if 4 < v < v c> \. 
If v > v Ct i then G(— 2+w; u) continues to decrease until it reaches zero and afterwards it becomes negative. 
This defines a second critical value v c ^ Eq + 2 such that = G(— 2 + v Cy 2', v c ,2) > G(—2 + v; v) for every 
v > f c ,2- This generates a non-degenerate eigenvalue X(v) in the interval (2, —2 + v). 

To summarize: if i> £ [0, u c ,i) there exists a unique non-degenerate eigenvalue in the interval (v+2, 00). 
If u C) i < w < u C! 2 there are no discrete eigenvalues, and if w c .2 < w then we again have a unique non- 
degenerate eigenvalue in the interval (2, —2 + v). 

There are no embedded eigenvalues in the set (—2, 2) U (—2 + v, 2 + v), for all v > 0; the explanation 



is that the imaginary part of G(x + i0+; v) is not zero if £ € (—2, 2) U (—2 + v, 2 + v), see (2.7) and (2.5 1 . 
The only remaining situation where eigenvalues could exist is at thresholds, i.e. when v is either v Ct i or 
v c ^2- But we will show in the next subsection that this is not the case. 

2.2 Behavior near thresholds 

Let us assume that v c .\ — v > is very small, which means that the eigenvalue \(v) £ (v + 2, 00) still 



exists but < X(v) - (v + 2) << 1. From (2.7) and (2.5) we obtain: 

F{t,v) := G(t 2 +v + 2;v), G(x;v) = F(y/x - v ~ 2, v), (2.11) 
F{t,v) = E -v-2-t 2 + ^T 2 (2 + < 2 - ty/A + t 2 + t 2 + v + 2 
- v/(w + t 2 )(i; + 4 + t 2 )), 

which admits a smooth extension near the point (0,u Cj i), with F(0,v cA ) = and d t F(0,v cA ) = -t 2 ^ 0. 
Then the implicit function theorem gives us a smooth map t(v) defined in a neighborhood of w c ,i where 
F(t(v), v) — 0. Since d v F(0, u c l ) = — 1 + 0(t 2 ) 7^ if r is small enough, we have t(v) ~ v c l — v > 0. It 
follows that 

X(v) =t 2 (v) + v + 2, (2.12) 
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and X(v) — (v + 2) ~ (u c l — v) 2 near the threshold. 
Moreover, since G(x; v) = F(y/x — v — 2, v) we have: 



d x G(x;v)~- - 0<x-(v + 2)«l, 

2y/x- (v + 2) 



(2-13) 



and this estimate holds for v in a small neighborhood of v c i. In particular, for v — v c \ and using that 
G(v c< \ + 2; v Ct i) = we have after integration: 



G(x;v c ,i) ~ -t 2 \Jx- {vci +2), 0<s-(» e ,i + 2)«l. (2.14) 
Proposition 2.1. If v is either v Ci i or v c ^., then neither v Ct i + 2 nor w c .2 — 2 are eigenvalues of h. 

Proof. We only give the proof for v Ct x; the other case is similar. We will prove that for every basis vector 
ip e {|m±) : m > 0} U {\S}} we have 



lim (x-v cl -2){ip\R(x)i>) = 0. 

2J\i"c,1+2 



(2.15) 



Coupling this with the fact that we always have 



lim (x - v c ,x - 2){i)\R(x)P ac i>) = 0, 



;P, then 



it shows that if there exists some orthogonal projection P such that R(x) = i?(x)P ac + , , 

necessarily (tp\Pip) — \\Pip\\ 2 — for all basis elements, hence P = 0. Let us prove (2.151 for ip = |m_). 
From (2.4) we have: 



(m_ |i?(x)m_) 

= (m_ |r_ (x — u C! i)rn_) + 



G(x;w Cj i) 



|(m_Mx-« Cil )0_)| 5 



(2.16) 



and from (2.6 I : 



; \ \ n , l-(Ci 2 (^-^.i)) m+1 
(m_|i?(x)m_) = -Ci(a; - w c ,i) — > — 

X (,j (X W c> iJ 



G(x; v c 



-\(i{x ~ v c< x)\ 



2m+2 



(2.17) 



Since Q (x — w c ,i) goes to 1 when x converges to w Ci i + 2, it follow that the only singular behavior comes 
from \2.\A\ , thus: 

(x — v Ci i ~ 2){m-\R(x)m-) ~ ^/x — u c ,i — 2 — > 0. 



The other cases are similar and we do not treat them. 



□ 



Proposition 2.2. The one dimensional projection Pd(v) associated with X(v) is differentiable on (0, u c ,i)U 
(v c 2 ,oo), and for every compact K C R i/iere exists a constant Ck > swc/i that we have the (optimal) 
bound: 

\\p»\\ < % - y «Gi^n[(0,« Cll )u( Uc , 2 ,co)]. 

Proof. We only concentrate on the case in which < v c i — v « 1, that is just before the eigenvalue 
X(v) € (w + 2, oo ) hits the first threshold and disappears. The eigenvalue X(v) — t 2 (v) + v + 2 is smooth 
as a function of v because t(v) is, see the discussion preceding (2.12). Using the Riesz formula for the 
projection Pd(y) associated with X(v), we see from (2.4) that it will consist of a finite sum of rank one 
operators, like for example 



Pa{v) 



(d x G)(X{v)-v) 



r_(2 + i>))|0_)(0_M2 + t>)). 



The above operator turns out to be the most singular when v lies near the critical values. From (2.11) 
we see that (d x G)(X(v); v) = ^/ff hence: 



2tH(v) 



(d t F)(t(v),v) 



r_(2 + t 2 ( V ))|0_)(0_|r_(2 + i 2 ( V )). 



(2.18) 
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If fc > 1 and \x\ > 2 we have: 

\\r k ± (x)\0 ± )\\ = ^<0±|4*(*)0±> 



V(2fe - 1)! 



^ fe - 1 (0 ± |r ± (x)0 ± ) = 



V(2fc-1)! ' 



(2.19) 



From (2.5) we see that 



which implies: 



In particular: 



ll^(x)|o ± )|| 



(x 2 — 4)2fe-3/2 



(x 2 - 4) fc - 3 /4 ' 



0<|x|-2«l, fc>l, 



0<|ac|-2«l, fc>l. 



(2.20) 



|r*(2 + t 2 («))|0±)|| 



i(«)2fc-3/2 

7 , 2A . , /2 , < v c ,i - 1> « 1, fc > 1. 

(D c X — u) ' 



(2.21) 



If fc = 1 we obtain that 1 1 r-__ (2 + t 2 (v))\0±) \\ 2 ~ l/t(u), which shows that P a (v) is bounded near the 
threshold. Now by differentiating ( |2 . 18 1 with respect to v, and keeping in mind that t'(v) is bounded, we 
see that the singular behavior is given by: 



max{^) 2 ||r2(2 + t 2 («))|0 ± )||||r_(2 + < 2 («))|0 ± )||, 

1 



M2 + t 2 ( W ))|0 ± }|| 2 } 



which is of the type claimed by the proposition. 

Now we can analyze all the other terms given by the Riesz formula, and notice that they contain 
at most one resolvent r_(2 + t 2 (v)), which is the only object which produces singularities. Since P a (^) 
contains two such resolvents, it is the most singular object. Clearly, for v > v Ct 2 the singular object would 
be T-(— 2 — t 2 (v)) and P a (v) is the term containing two such resolvents. It turns out that near both 
critical potentials v c we have: 

ll-Pd(v) - P a (v)\\ = O(VF^I), (2.22) 
\\P' d {v)-P' a {v)\\=0(\v~v c \- 1/2 )- (2.23) 

We now prove that Pd(v) converges strongly to zero when v — ¥ v Cl i. First, if m > is fixed, then from 
(12.180 and P^l we have: 



(m_|P a (w)m_) = - 



2tH{v) 



C 2 1 m+2 (2 + t(v)'')^t{v). 



(d t F)(t(v),v) 



Second, using Pd(v) — Pd(v) and (2.22) we have: 



\\P d {v)m^f < ||P d («) - P a («)|| + (m-\P a {v)m-) ~ y/v cA - v. 

Up to an e/2 argument one can now show that Pd(v) converges strongly to zero. All this is consistent 
with the fact that there are no eigenvalues at thresholds when v = v c . The projections Pd(v) delocalise 
more and more and converge strongly to zero. Notice that this implies that if Pd(v)ip(y) — ip(v) then 
■0(f) converges weakly to zero. □ 



3 Proof of Theorem 1.1 



We only prove in detail (1.10), the other identity follows after a similar argument. 
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3.1 A general propagation estimate 

Let s_i < so < and consider a function v : (—00, sq] <— > R such that u(s) — ^(s_i) for s < s_i, 
and v is piecewise continuous, in particular continuous at sq. The values v(s-i) and z^(so) are n °t 
one of the two critical values i> Cj i and tj C 2- Denote with P ac j (respectively -Pd,i) the projection on the 
absolutely continuous (respectively discrete) subspace of h + v(s_i)H_, and with P ac { (respectively P<i,f) 
the projection on the absolutely continuous (respectively discrete) subspace of h + ^(s )II_. 
Remember that ho = hs + h^ and h = Hq + hx- The subspace of absolute continuity of ho is 

P ac (h L ) = o®n_®n + . 

Denote by 

Q;:=s lim e «(h+»'(«-i)n-) e -«(ho+K-- 1 )n_)p ac ^ J .\ ( 3A \ 
t— 00 

the 'incoming' wave operator between ho + ^(s_i)II_ and h + i^(s_i)n_. This operator exists and is 
unitary between Ran(P ac (/iL)) and Ran(P aC; f) IB] - In a similar way, define: 

ilf.= s lim e it(h+ v (s )n.) e -it(h +u(s )n.)pf hL Y (3.2) 
t— 00 

If t < w < so/77 we denote by U(t,w) the unitary solution to the evolution equation 

iU'(t,w) = h(r]t)U(t,w), h(s) := h + v(s)Il-, (3.3) 
U(w, w) = 1, —00 < t < w. 

Proposition 3.1. Let $ e W de a fixed unit vector. Using the notation Au :— i v(t)(1t. we have the 
following estimates: 

lim c/(£zl ^)*r) ie - ji ^P^e^ Al/n -P ac (/i L )* = Q/P^)*, (3.4) 



which is equivalent with: 



Iim||l7(^=i,^)n / P ac (fe £ )* 

(«-l-a ) , 



^e" 4 ^r— ^e^n-p^^)*!! ^0 (3.5) 
and 

lim||C/( S ^, S °)P ac , f vI/ 
v\fi V 77 

- ae~ lll ^P l2 ^ e l ^r n -p ac (/i L )mp ac , f *||= 0. (3.6) 



Proof. We start by proving (3.4). Denote by Uo(t, w) the explicit unitary solution of the equation 

iUo(t,w) = (h L +v{rjt)TL_)U (t,w), U (w,w) = l, -00 < t < w 
U (t,w) = e-^~ w )^ e ~^^> (s)dsn -. (3.7) 

Note that Uo(t,w) commutes with P ac (/iL)j and equals the identity operator when restricted to the 
subspace of the sample. 

Choose / € RanP ac (/i £ ) such that if is a generalized eigenfunction of /i £ , then the function /) 
is a smooth function of E with a compact support not containing ±2. First, (5, C/o(i, w)f) — 0. Second, 
from (3.7) and by partial integration with respect to E it follows that there exists C/,at > independent 
of 77 such that 

max :{\(j ± ,U (t,w)f)\} < C ftN (t - w)' 2 , Vt<w. (3.8) 

0<j<N 

Define 

W(i):=^(i,-)*^o(<,-)P«(^), i<-. (3.9) 
77 77 ry 
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We will compute the strong limit of W(t) when t —> — oo in two different ways. Up to an e/2 argument 
it is enough to prove the existence of a strong limit for a function F such that F) is a smooth function 
of E with a compact support not containing ±2. 

On one hand, if t < S—i/t) we have: 

W(t)F = U(t, ^yu (t, -)P ac (h L )F 

= U( — , ^)* e i(t-s-i/v)h(v(s-i)) 
^ T) ' T) 

■ e'^- s "/^ P ac (h L )e~ * ^ dsYI - F 

= Tj( ! L±^ f°")* e i(«-s-i/'7)h( ! /(s_i)) e -i(t- S - 1 /j7)(/Jo+^(s-i)n_)p ac ^ i N 
7] ' T) 



-Ai>n_ 



F. 



(3.10) 



• e i e i 

where in the second line we used the group property of U and the fact that v is constant at the left of 
S-\. The variable t only appears in the middle, and using (3.2) we obtain: 

,s_i s 



lim W(t)F=U( — ,—) Qi 

t— y— oo *? ^? 



P ac (h L )F. 



(3.11) 



On the other hand, by differentiating with respect to t the formula (3.9) defining W(t)F and then 
integrating back we obtain: 



r-t 

W{t)F = P ac {h L )F + i / dTU(r,^Yh T Uo(T,^-)P ac (h L )F 
= P ac {h L )F 



+ i 



drU(r + S -°, ^Yh T Uo(r + ^, ^)P ac (h L )F, 
r\ rj fj Tj 



so S x 



(3.12) 



or using (3.11 1: 



lim W(t)F = U(?=± S ^Ya i e- ll ^i fRh -e^ A ' Jn -F = P ac (h L )F 

t- — y — oo 7] Tj 



+ i 



drU(r + ^)*h T U Q {r + ^)P ac (h L )F 
rj rj rj 77 



For fixed r, we have: 



= l + i T dt'Ult' + ^)*(i/faf + s ) - K S o))e- lt '^ (so)) , 
Jo V V 



which leads to the operator norm estimate: 



(3.13) 



|tf( T+ ?E ^)V <Th < v < ao »-l||<T sup |i/(ar) - e/(« ) 
77 77 



(3.14) 



or 



|C/(r+ ^, ^) - e -^(KV»|| < T sup | v ( a )-i/(s )|, 

?7 Tf i)T+sa<x<s 



(3.15) 



where the right hand side goes to zero with 77. The same type of norm estimate as (3.15) holds when we 
replace U with Uq and h with hi,. Being a norm estimate, it also holds for the adjoints. 
Using (3.8), we obtain the norm estimate: 



\h T U (r + ^, -)P ac (h L )F\\ < C(t) 
1 77 77 " 



(3.16) 
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Now using the left continuity of v at s , the L 1 upper bound in ( 3.16 ) and the pointwise norm convergence 
of (3.15), the Lebesgue dominated convergence theorem gives us (3.4): 



= P,c(h L )F 

-oo 



P ac (h L )F 



= n f p ac {h L )F, 

where the last equality comes from the Dyson equation satisfied by flf. Then since U(^-, ^) is unitary, 
we have: 



Ui^^SlfP^h^F-Cke 
v rj T) 



- hL e$ AvTl -P aK (h L )F = o{l) 



which is exactly (3.5). Finally, if we use in (3.5) a vector F — f2£P aC) f\& where is some arbitrary unit 
vector, then Q,fP ac (h,L)F — P aC: i^ and (3.6) follows. 



□ 



3.2 Proof of (1.10) 



Proof. Note that it is enough to consider rank one observables, thus we will assume without loss of 
generality that A = |^)(^|. Let us remember the simplifying notation: 

Pac(0) := P ac (h(v(0))), f a c(-l) := Pac(M«(-l))) = P ^(h), 

P d (0) := P d (h(v(0))), P d (-1) := P d (h(v{-1))) = P d (h). (3.17) 

We can write: 

U(-r 1 -\0)*f eq (h)U(-r l - 1 ,0)A 

= u(- v - 1 1 oyf cq (h)p d (-i)u(- v -\o)p d (o)A 

+ U(-r)- 1 ,0)*f eq (h)P ac (-l)U(-r ] - 1 ,0)P a . c (0)A 
+ U(-r / - 1 ,0)*f cq (h)P d {-l)U(-rr\0)P,c(0)A 

+ U(-r ! -\0)*f cq (h)P lic (-l)U(-ri- 1 ,0)P d (0)A. (3.18) 
We will treat these four terms separately. 

3.2.1 The discrete-discrete term 

From now on the notation T\ ~ T 2 will mean that the difference T\ — T 2 goes to zero with 77 in the appro- 
priate topology. Under our assumptions on v, the discrete eigenvalues of the instantaneous Hamiltonian 
h(s) = h(v(s)) remain at a positive distance from the continuous spectrum if —1 < s < s c and s' c < s < 0. 
Moreover, we assumed that v"(s) is continuous and bounded on the open intervals. Then using the usual 
adiabatic theorem for eigenstates, the continuity and the group properties of U, we obtain that in the 
operator norm topology we have: 

P d (-l)Z7(-l/r?, s c /r?) ~ U(-l/r], s c /^)Pd(s c ,i - 0), 
U(s'J V , 0)P d (0) ~ P d (s' c + 0)U(s'Jri, 0), 
P d (-l)U(-r,-\0)P d (0)~ 

U(-l/r], s c /r))P d (s c ~ 0)U(s c /ri, s' c /ri)P d { s ' c + 0)U(s' c /r}, 0). (3.19) 

We see that in the middle of the last line we get the rank one operator 

P d (s c - 0)U(s c /r), s'Jrj)P d {s' c + 0). 

We will now show that 

lim(^(s c - 0) I U(s c /ri, s'Jr,)^ + 0)) = 0. 



In order to do this, let us go back to Proposition 3.1 and identify s_i = s c , s = s' c and v with v restricted 



to [s c ,s' c ]. We see that the potential v is continuous from the left at s' c , and moreover, the spectrum of 
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the instantaneous Hamiltonian h(s) is purely absolutely continuous if s € [s c ,s' c ]. Remember that the 
wave operators fl(s' c — 0) and 0(s c + 0) defined in (1.7) map Ran(P ac (/ii)) onto %. Comparing with 
(3.6), we can identify 0/ with Q(s' c — 0) and Oj with 0(s c + 0), while RanP aCi f = %. Thus: 

(ip{s c -0)\U (s c /<n, s'Jrj)ip{s' c + 0)) 

~ (0(s c + 0)>(s c - 0) | e- lS -^ hL e^ n -n(s' c - Q)*>(s' c + 0)). 

The rj dependence on the right hand side is explicit, and this term converges to zero with 77. It means 
that there is no contribution to the final steady state from the purely discrete part. 

3.2.2 The continuous-continuous term 



The next term in (3.18) contains both the initial and final projections on the instantaneous absolutely 
continuous subspaces. Note the important thing that P ac {— 1) = P ac (M commutes with f cq (h). 

We can again apply Proposition 3.1 where now s^i — —1, s = 0, O,- = 0(— 1), 0* = O(0) and 
P ac j = Pac(0). Using piS]) we get: 



^(-»r 1 ,o)*/c q (M^ac(-i)c/(- f ]- 1 ,o)p ac (o)* 

~ U(-r)-\ 0)*/eqWPac(-l)O(-l)e^ 

P ac (fc L )O(0)*P ac (0)*. 



e 1 



(3.20) 



The important thing is that 0(— 1) intertwines between ft, and h^. Thus the right hand side of the above 
equation is equal to: 

c/(-7r 1 ,o)*Pac(-i)o(-i) e ^- e i ^ n -/ cq (ft L )p ac (ft L )o(o)*p ac (o)*. 



(3.21) 



(3.22) 



But now we can apply ( |3.4[ ) and obtain the result: 

^(-^ 1 ,0)*/cq(M^ac(-l)C/(-^ 1 ,0)P ac (0)* 

-O(0)/ cq (ft L )P ac (ft L )O(0)*P ac (0)*. 
Thus the first nonzero contribution to (A) in ( |1.8| ) is: 

(O(0)*P ac (0)* I / eq (/lL)Pac(MO(0)*P ac (0)*). 

3.2.3 The mixed terms 

The third contribution in ( |3.18[ ) is: 

U(-v~\ 0rU q (h)P d (-l)U(- V -\ 0)P ac (0)*. 

A direct application of ( |3.6[ ) with s_i = —1 and so = shows that the vector U(— r; -1 , 0)P ac (0)4' lies 
almost completely in RanP ac (— 1), thus when projected on RanPd(— 1) it will converge to zero with 77. 
Thus this term will also disappear. 

The fourth contribution is probably the most interesting one. Using the decomposition 1 = P ac (0) + 
Pd(0) to the left, we can write: 



U{-v~\ 0)* /c q (h)P ac (-l)U(-f]-\ 0)P d (0)A 

= P ac (0)U(-i 1 -\0)*f cq (h)P llc (-l)U(-7 1 -\0)P d (0)A 
+ P d (0)U{-7 1 - 1 ,0)*f cq {h)P ac (-l)U(- V -\0)P d {0)A. 



(3.23) 



We show that the trace of the first term will converge to zero. After taking the trace with A = | , 3 / )( V E'| 
we obtain: 



(vE'|P ac (0)P(-7 ? - 1 ,0)*/cqWPac(-l)C/(-^ 1 ,0)P d (0)*) 

= ((7(-ry- 1 ,0)*/cqWPac(-l)C/(-^ 1 ,0)P ac (0)*|P d (0)*). 



(3.24) 



But the 'bra' vector in the second line is the same as the one in (3.21 ), which we know that asymptotically 
enters the range of O(0), thus becomes orthogonal to the 'kef vector. 
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Now let us treat the last contribution. Because -Pd(O) is one dimensional, the trace of this last 
contribution is: 



(^(0)|^(0)) 



Tr{P d (0)U(-r 1 -\oyf eq (h)P ac (-l)U(-r,- 1 ,Q)P d (0)}. 



(3.25) 



Using the group property U(—rj 1 ,0) = U(—i] 1 ,s' c /rj)U(s' c /r],0), then applying the adiabatic theo- 
rem which says that in operator norm: 

Uis'Jr,, 0)P d (0) ~ P d (s' c + 0)U(s'Jri, 0), 
P A {Q)U{s' c h, 0)* ~ Uis'Jr), 0yP d (s' c + 0), 

and finally the trace cyclicity, we obtain: 

Tr{P d (0)U(-r 1 - 1 ,0)*f cq (h)P llc (-l)U(- V -\0)Pd(0)} 



(<M4 + o)\u(-v~\ s ' c hTU{h)PU-^)u{-r\ s'JvMs'c + o)>. 



(3.26) 



Now we can reason as in the continuous-continuous case, where s_i = —1, sq = s' c , and Q{s' c — 0) plays 
the role of fi/, which ends the proof of (1.10). 

Regarding (1.9), the only difference is that the discrete-discrete term will now contribute because we 
can apply the usual adiabatic theorem on the whole evolution interval, while the mixed terms disappear. 
The continuous-continuous term is identical with the previous one. □ 



4 Proof of Proposition 1.3 



First of all, the result stated in (i) is a trivial consequence of the theorem. We will only prove (ii) in 
detail since the proof of (iii) is very much similar. 

Let us now focus on the case described in Figure [2j The family h(v) is norm continuous in v and 
its spectrum is Lipschitz continuous in v. The spectrum of h(v c ,i) is purely absolutely continuous and 
consists of two well isolated bands: [—2, 2] U [—2 + v Ci i, 2 + u c ,i]. We can find a positively oriented circle 
T which completely includes the interval [— 3 + v c .i, 3 + u c ,i] while [—2,2] lies outside. Moreover, if S > 
is small enough and \v — t> Cj i| < Sq then the spectrum of h(v) will be at a distance proportional with Sq 
from the spectrum of h{v Ct \) and: 

sup sup \\{h{ v ) ~ < 00 • (4-1) 

zer \v-v cA \<5 a 

The above bound also holds true if h{v) is replaced with h^v) = + vll_. 
Using the second resolvent identity, we have that 

U_{(h(v) - z)- 1 - (h L (v) - z)- 1 } = -Il_(h L (v) - z^hrMv) - z)~\ 
for zeT. We now argue that if = ^2 m>M |m_)(m_| is IT_ without its first M sites, then: 

lim {sup sup \\n {M) (h L (v)-z)- 1 h T (h(v)-z)- 1 \\} = 0. (4.2) 

M->oo zer |t;_ Vcil |<5 

This is due to the fact that is localized while, uniformly in z E T, the two resolvents have exponentially 
localized kernels near the diagonal. 

By integrating (/il(v) — z)^ 1 over T we obtain n_. Thus the Riesz projection P(v) corresponding to 
the spectrum of h(v) contained inside T obeys the estimate: 

lim { sup ||n (A/) { J p(w)-n_}||} = o. (4.3) 

M->oo \ V - Vcl \<s 

Now fix some e > 0. We can find M — M e such that: 

sup ||ny e) P(^)-nL Me) || < -„/„ . (4.4) 



\v-v cA \<S 



6||/ cq ||c 
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In order to shorten notation, let us write ifi c g instead of ip(s' c + 0). This is the eigenvector corresponding 
to the bias v = v c> \ — 6, and the eigenvalue X(v) is only barely larger than v + 2, their difference being 
proportional to <5 2 , see (2.12). 

From (2.22 1, (2.181 and the fact that ip c< s converges weakly to zero as 5 — > (see the end of Section 



2.2 1 we see that ip Ct g becomes more on more delocalized, far away on the left lead, when 6 tends to zero. 

(4.5) 



In other words, for every fixed M > we have: 



iim||^ c , 5 -nL M) ^MI =o. 

Thus there exists S £ small enough such that for every 5 < 5 e we have: 

T (M«) , „ - 



c,S 



it 



6||/< 



(4.6) 



cq 1 1 oo 



Again in order to shorten notation, let us write Q Ct $ instead of fl(s' c — 0). Going back to ( |1.11[ ) and using 
(4.6 1 we obtain: 

\{lpc,5 I ^c,i/eq(/li)-Pac(/lL)^c5' ( /'c,5) 



< 



which together with (4.4) and (4.6) gives: 



\(ipc,s I ^c,sf C q(h L )P ac (h L )n* c S Tp Ct s} 

- {ll> c ,s\P{v)n c , S f cci {h L )P &c {hLW c ^c,8)\ < 



(4.7) 



for every \v — v c ,i\ < <^o and 5 < 6 e . 

The idea is to take v — v c .\ + 5 in the above estimate, corresponding to v(s' c — 0) for which the 
Hamiltonian has purely absolutely continuous spectrum. Before that, let us notice a few identities. Using 
the intertwining property of Q Ci $ we have that P(v c> \ + 5)Q C $ = i7 Ci 5iI_ where II_ is nothing but the 
spectral projection of the decoupled operator /io + (^c,i+^)n_ corresponding to the absolutely continuous 
spectrum contained in T. Defining f e q,s(x) ■= feq(% ~ v ci ~ 5) we have the identity: 

Ti-f cq (h L ) = n_/eq, a (/»£ + (vd + <y)n_), 

or putting everything together: 

p( VcA + S)n c , s f cq (h L )p ac (h L )ni s 



= p{v c ,i + 5)n c , s f cq j(h L + (u C)1 + 5)n_)fi* ;5 

= P{v c ,i+ 5)f cq , s (h(v cA +5)), 



(4.8) 



where the second identity is again implied by the intertwining properties of il c ,5i remember that h(v c ,i+5) 
has purely absolutely continuous spectrum. 

Thus putting v = u C) i + 5 in (4.7) and using the identity (4.8), we obtain that for every 6 < S e we 
have: 



\{ipc,5 I Qcsfc^h^Pacih^n^gip^s) 

- {i>c,s\P(v c ,i +5)f cq ,s{h{v cA +S))^g)\ < £ -. 



(4.9) 



Using again (4.6) and (4.4) we can replace P(v Cl i + S) in (4.9) with the identity operator, with a cost of 
an error of at most e/3. Thus for every S < 5 e we have: 



I (ipc,S I Qc,8f cq (h L )P ac (h L )n* s ip c j) 

- (i>c,S I feq,s(K V c,l + S))lp c j) \ < 



he 



(4.10) 



Since the family h(v) is norm continuous and f eq is supposed to be continuous, we may find a 5' e < S e 
such that 

\\f cq AHvc,i + 6)) - U q ,s{h{v c<1 - 8))\\ < e/100, S < 5' e . 
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But tjj Ct g is the eigenvector of h(v c _i — S) corresponding to the eigenvalue A(v Cj i — S). Thus: 

\(lpc,S \^sfcq(hL)P a c(hL)tt* c S 1p c ,s) 

- /eq(A(t> c ,l -6)- « c ,i - *)| < y , V<5 < 

But A(w C: i — 5) — v c ,i — S converges to 2 when 5 goes to zero, thus there exists S'^ < S' £ such that for every 
S < S s we have: 

\(fl>c,6 I ttc,6fe q (hL)P a c(hLW c ,si>c,6) - /eq(2)| < £. 

The proof is over. □ 
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